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Abstract. We describe the electrically enhanced growth of needle crystals from 
. the vapor phase, for which there exists a morphological instability above a 
^ threshold applied potential. Our improved theoretical treatment of this phe- 

Q . nomenon shows that the instability is present in both polar and non-polar sys- 

o: 

terns, and we provide an extension of solvability theory to include electrical 

>- . effects. We present extensive experimental data for ice needle growth above the 

^ . electrical threshold, where at T = —5 C high-velocity shape-preserving growth is 

^ . observed. These data indicate that the needle tip assumes an effective radius R* 
^ \ 

O . which is nearly independent of both supersaturation and the applied potential. 
> I 

g . The small scale of R* and its response to chemical additives suggest that the 

. needle growth rate is being limited primarily by structural instabilities, possibly 

O . related to surface melting. We also demonstrate experimentally that non-polar 
systems exhibit this same electrically induced morphological instability. 
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The formation of stable spatial patterns is a fundamental problem in the study of nonlin- 
ear nonequilibrium systems and controlling pattern formation has generated considerable 
recent interest in light of a host of possible technological applications. A now-standard ex- 
ample of a condensed-matter pattern-forming system is the diffusion-limited growth of free 
crystalline dendrites, which are nearly ubiquitous products of rapid solidification, from ei- 
ther liquid or vapor precursors. While the diffusion equation alone is sufficient to define a 
relationship between the dendrite tip velocity and tip radius, typically surface tension (which 
changes the surface boundary conditions via the Gibbs-Thomson mechanism) must be in- 
cluded in order to select a unique needle-like solution. Microscopic solvability theory has 
succeeded in furnishing a mathematically consistent and dynamically stable solution to this 
problem for simple 2D and 3D dendrite growth (ignoring, for example, surface kinetic and 
surface transport effects) Instabilities and noise amplification leading to sidebranch 

generation have also been well studied. 

In this Letter we examine shape-preserving needle growth from the vapor phase in the 
presence of an applied electrical potential. We previously described how high electric fields 
and field gradients near the needle tip can enhanced its growth, and subsequently drive in a 
new kind of growth instability for this system [^,?,?,?]. The growth behavior for small ap- 
plied potentials is quite well described by an extension of solvability theory which we describe 
below. Above a threshold potential, however, the needle growth cannot be stabilized by the 
Gibbs-Thomson mechanism, and so the phenomenon lies outside the realm of solvability 
theory. We performed a series of experiments to examine needle growth in this new electric 
regime, where we can still observe shape-preserving growth, but with much increased tip ve- 
locities.. From these measurements we are able to infer the physical mechanisms responsible 
for controlling the enhanced needle growth rates. 

We consider the case of diffusion-limited growth of a free crystalline dendrite, which is 
well known to exhibit a constant- velocity, needle-like solution with a parabolic tip geometry 
0]. We will assume growth from the vapor phase, with low solute concentration in a non- 
polar solvent gas. When an electrical potential is applied to the growing crystal, the normal 
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diffusion equation is replaced by the Smoluchowski equation [|,?] 

|^ = D^.(^c + c^$) (1) 

where c is the solute concentration, D is the diffusion constant, and the external force felt 
by the solute molecules is described as the gradient of the potential 

where a is the molecular polarizability, and the electric field is the gradient of the electrical 
potential ^ = — "^^9 ||^. Ignoring interface kinetics, the continuity equation at the interface 
yields the normal component of the surface growth rate Vn = in ■ if surf) as 

Vn = ■ {"^C + c'^^^urf (3) 

C-solid 

where CgoUd is the solid number density and the right-hand side is evaluated at the solidifica- 
tion front P). For our experiments the diffusion length ^£, = 2D/vn is very large compared 
to other length scales in the problem, so we can assume the slow-growth limit where dc/dt 
can be taken equal to zero. 

The boundary conditions in this problem are defined by the equilibrium vapor pressure 
above the crystal surface, which depends on the applied electrical potential, along with 
the vapor pressure at infinity, Cqo- For a spherical crystal with radius R a perturbation 
calculation yields the vapor pressure 

+ f) (4) 

where Csat is the equilibrium vapor pressure above a fiat interface, = 2[3 / CsoUdkT, (3 is the 
surface tension, i?^^ = eo(pl/2csoHdkT, R^^i = aip^/kT, and (po is the applied potential 0. 
The second term in the above expression is the usual Gibbs-Thomson effect 0, the third 
term arises from the electrostatic energy of a charged sphere, and the last term is present 
for polar molecules, since we are assuming zero field inside the crystal. 



3 



With these modifications to the diffusion equation and surface boundary conditions, we 
then obtain a modification of the Ivantsov equation relating the tip velocity, and tip 
radius, R, of a parabolic needle crystal P,?] 



2-D Csat 

R\og{r]oo/R) Csoiid 



elec 



R i?2 



(5) 



with 



'''' CsomkT\og\ii^/R) 

where Ai = (coo — Csat)/csat, C = 1 + 2'yaCsoiid^i/^o, and 7 is a dimensionless geometrical 
constant. From the spherical growth problem, solved in the limit of small Reiec using a 
perturbation approach, we expect 7 ~ 0.2. For (fQ = this becomes the slow-growth limit 
of the Ivantsov solution if we identify rjoo ~ ^d- This result supersedes previous treatments 
0,?,?] in which a cruder approximation was used, which effectively assumed 7 = 0. For ice 
we find C ~ H-4.7Ai, and Reiec ~ 40(v9o/1000 V) nm 0; thus in our experiments described 
below the polarization term is larger than the electrostatic term. For unpolarized molecules 
C = l. 

In the absence of an applied potential solvability theory tells us that the small capillarity 
term, (Iq/R, stabilizes the needle growth; the decrease in growth velocity for small R brought 
about by the Gibbs-Thomson effect opposes the MuUins-Sekerka instability |jlO| . The applied 



potential acts to increase the tip velocity for smaller R, thus tending to destabilize the 
growth. The electrical perturbation can be incorporated quite simply into solvability theory 
by defining a modified capillarity, (Iq = do — CR^i^jR, which then yields the new solvability 
relation 

,„ = ^R£^(l-C§!§) (7) 

vR^ Csoiid V doR J 

where (Tq is the solvability parameter in the absence of an applied potential. Combining this 
with the Ivantsov relation then yields a quadratic equation for the tip radius, R^ — RqR + 
CRoRli^Jdo = 0, where Rq is the tip radius in the absence of an applied potential. This 



quadratic equation has no real roots for Reiec > Reiec,thresh = (doRo/iC) ' , implying that at 
above a threshold potential the Gibbs-Thomson effect can no longer stabilize needle growth. 

Ice is a convenient experimental system in which to observe this phenomenon, and we 
previously described the behavior of ice needle growth for applied potentials at or below 
threshold @,?,?]. The existence of a threshold potential was established, and the tip velocity 
as a function of (^Jq below threshold is described by the above theory (although there remain 
considerable theoretical uncertainties associated with the growth of faceted crystals). A 
principal aim of the experiments described here is to quantitatively examine the growth 
behavior above threshold, beyond the realm of solvability theory. 

Our experiments were performed in air at atmospheric pressure, using a vertical diffusion 
chamber in which the supersaturation could be changed by adjusting the temperature profile 
within the chamber [Q. Crystals were grown on a thin wire substrate, to which an electrical 
potential could be applied. Supersaturation was determined using a combination of diffusion 
modeling and frost-point measurements, and we estimate an overall scaling uncertainty in 
our reported Ai of roughly ±30 percent. 

Crystals grown at T = — 15 C exhibited a classic dendritic morphology with growth 
along the a-axis, and the tip velocities were observed to increase approximately linearly 
with Ai as shown in Figure 1(a). Above threshold at this temperature we nearly always 
observed the tip-splitting phenomenon described in [^]. Crystals grown at T = —5 C were 
needle-like, and at these high supersaturation levels the needle axis was typically displaced 
from the c-axis, with larger angular displacements toward the a-axis observed for larger Ai. 
Fits to these data using the Ivantsov relation yield tip radii of -Rq ~ 1-2 /um and 1.5 /xm 
for T = —15 C and —5 C, respectively, independent of Ai. In both cases the data do not 
follow the trend f ~ Af expected from solvability theory, at least over the limited range in 
Ai shown here. This is probably because the crystals exhibit some faceting, implying that 
growth is limited by surface kinetics in addition to diffusion. Solvability theory does not 
extend to the faceted case, but we expect that the additional complication of facets will not 
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affect our main results to a very large degree. 

Above threshold the tip velocities of T = —5 C needles increased substantially, as is 
shown in Figure 1(c) and Figure 2. The threshold potential was found to be approximately 
1000 volts, independent of Ai, implying Reiec,thresh ~ 40 nm. This is comparable to the 
expected (cio^o/4C)^/^ ^ 30(1 + 4.7Ai)-i/2 

nm (using = 2 nm and the Rq inferred from 
our measurements). The quantitative agreement is satisfactory, but again the observations 
do not show the supersaturation dependence suggested by our simple solvability theory. 
From the supersaturation dependence at (po = 2000 volts we find tip radii of R* ~ 360 nm 
using the Ivantsov relation, again independent of Ai. These growth rates were observed to 
be the same (to ±15 percent) for positive and negative applied potentials. 

Adding trace quantities of certain chemical additives to the air in our diffusion chamber 
changed the electric needle growth dramatically, as shown in Figure 1(c). Under certain 
conditions the electric needles then grew along the c-axis, with velocities approximately four 
times faster than without additives, whereas there was no perceptible change to the normal 
crystal growth |TT|. This chemically and electrically mediated growth as a function of (fo 
was qualitatively similar to the behavior shown in Figure 2 scaled by a factor of four, but 
the needle-to-needle scatter in the data was considerably greater. These data imply R* ^ 90 
nm at <^o = 2000 volts. 

With these data in hand we can then examine possible physical mechanisms responsi- 
ble for limiting the Mullins-Sekerka instability and producing shape-preserving growth at 
a constant tip velocity. Latent heat generated at the growing tip is a possible stabilizing 
mechanism, but it appears that the resulting thermal effects are not very important. As- 
suming heat generated by condensation is conducted along the needle and into the solvent 

— 1 /2 

gas, then the temperature rise at the tip is approximately AT ^ {t^soHdi^soivent) LpvR, 
where L is the latent heat and p is the solid density 0. Using the Ivantsov relation for vR, 
this becomes AT ^ O.O6A1 °K. This temperature increase is small, and does not exhibit the 
strong R dependence needed to stabilize the tip growth. 

Field emission can be examined using the Fowler- Nor dheim relation, from which we esti- 
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mate a current (for y^o < 0) of / ^ 14mA-((/?o/1000 V)2exp [-310(i?/l ;um)(1000 V/ |^o|)] • 
Assuming a power P ~ lipo is deposited at the needle tip (likely an overestimate), and 
assuming the same simple conduction model as above, we estimate that tip heating be- 
comes significant for / > 10^^^ A. Ice surface conductance is high enough |[T^ that such low 
currents produce a negligible voltage drop along a needle. 

We measured the current emitted from growing needle crystals, for which there is a 
considerable dependence on the sign of the applied potential. For ipo = 2000 volts we 
find / < 1 pA in clean air and / 1 pA for chemically mediated growth (with considerable 
scatter in the latter measurements). For (po = —2000 volts we measure 1^1 — 3 pA without 
additives and J ^ 10 pA with additives (again with considerable scatter). When / > 10 
pA we observe that the needle growth slows substantially, and the tip velocities become 
quite variable from needle to needle. We conclude that field emission plays a significant role 
when (fo < —1000 V and the crystals are grown in the presence of chemical additives, which 
results in sizable currents. In all other circumstances it appears field emission does not limit 
the needle growth. 

Having eliminated thermal and field emission mechanisms, we propose that the electric 
needle growth is stabilized by the diminished structural integrity of the material when formed 
into an extremely sharp tip. That is, the needle tip radius becomes so small that its solid 
structure is unstable and subject to deformation driven by surface tension over timescales 
comparable to the crystal growth time ~ R/v. There is considerable uncertainty in our 
understanding of the structural properties of ice, however, particularly for such extremely 
small crystals, which makes detailed calculations difficult [0. Surface melting certainly 



sets a lower bound on tip radius, since surface melting is known to produce a quasiliquid 
layer at these temperatures, which likely has very little resistance to shear stresses. The ice 



quasiliquid layer has a thickness of roughly 20 nm at T = — 5C [jT3[, which is not insignificant 
compared to the tip radii inferred from our measurements. 

One feature of our data that is readily explained by this model is that the measured tip 
radius R* is essentially independent of Ai and ipo in the electric growth regime (at least 



above a Ai-dependent transition region; see Figure 2). Above the threshold potential the 
tip radius decreases and the tip velocity increases to the point that structural deformations 
become a dominating influence. This must occur at some material-determined size scale, 
so we would expect that the new tip radius would depend mainly on material properties, 
and would not depend strongly on supersaturation or the applied potential, as is observed. 
The change in R* in the presence of additives could be explained by the strongly anisotropic 
structural properties of ice, since the most rapid growth observed is for c-axis needles. We 
see that this stabilization mechanism is a distinct, more direct, manifestation of surface 
tension. While normal needle growth is stabilized by the increased vapor pressure over 
surfaces with large curvature — the Gibbs-Thomson mechanism — electric needle growth is 
stabilized by structural deformations brought about by surface tension in the sharp tip. 

In a separate experiment, we also examined the electric growth phenomenon for a non- 
polar molecular solid by examining the growth of iodine crystals from the vapor. At a 
supersaturation of Ai ^ 0.1 we obtained the results shown in Figure 3. Rapid needle 
growth above a threshold potential of (fo ~ 1000 volts was clearly observed, although there 
was considerable scatter in the data. These data provide qualitative support of the above 
theory, and clearly verify that this morphological instability is present in non-polar systems. 

In summary, we have examined the detailed physics of electrically enhanced needle crystal 
growth from the vapor phase. We flnd that an extension of solvability theory explains the 
existence of a threshold applied potential, ipthresh, and adequately reproduces the magnitude 
of f thresh as wcU as the needle growth behavior for (po < (pthresh- For (po > ipthresh we 
observed high-velocity growth of ice and iodine needle crystals, demonstrating that this 
new growth instability is indeed present in both polar and non-polar systems. Extensive 
measurements with ice suggest that electric needle growth in this system is usually not 
limited by the effects of latent heat deposition or fleld emission, but appears to be limited 
by structural deformation of the needle tip. Additional experiments with other materials 
would further elucidate this process. In particular, since ice is a relatively soft material W2 



these investigations suggest that harder materials could yield much sharper tips. Related 



experiments with refractory metals have yielded electric needle crystals with nanometer- 
scale tips [|n|, and we believe that additional research into electrically enhanced growth 
may lead to useful applications. 
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